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We consider the statistical properties of vacua and inflationary trajectories associated with a 
random multifield potential. Our underlying motivation is the string landscape, but our calculations 
apply to general potentials. Using random matrix theory, we analyze the Hessian matrices associated 
with the extrema of this potential. These potentials generically have a vast number of extrema. If the 
cross-couplings (off-diagonal terms) are of the same order as the self-couplings (diagonal terms) we 
show that essentially all extrema are saddles, and the number of minima is effectively zero. Avoiding 
this requires the same separation of scales needed to ensure that Newton’s constant is stable against 
radiative corrections in a string landscape. Using the central limit theorem we find that even if 
the number of extrema is enormous, the typical distance between extrema is still substantial - with 
challenging implications for inflationary models that depend on the existence of a complicated path 
inside the landscape. 


I. INTRODUCTION 

We study the statistical properties of potentials with 
many degrees of freedom, and the inflationary trajec¬ 
tories that exist inside them. A large part of our in¬ 
terest in this problem follows from the string landscape 
[1-3]. However, we focus on the statistical properties 
of arbitrary multifield potentials, without making direct 
appeals to stringy physics or supersymmetry. Conse¬ 
quently, there is no guarantee that the string landscape 
will conform to the conclusions we reach. Rather, this 
paper determines what we can expect solely on the basis 
of the string landscape being a large, multidimensional 
potential.^ In doing so we have two aims. Firstly, to the 
extent that the properties of the string landscape are a 
function solely of its dimensionality and other simple pa¬ 
rameters, the analysis presented here elucidates its prop¬ 
erties. Secondly, if the actual string landscape does not 
match the conclusions we reach here, we will learn that 
its properties follow from the structure of string theory 
itself, and not just large-A^ statistics. 

A program for computing the number of vacua in the 
string landscape has been implemented [5-8]. However, 
in keeping with our interest in what follows from the 
large dimensionality of the landscape, our starting point 
will be an elegant argument championed by Susskind 
[4], which explains why we expect a landscape to pos¬ 
sess a huge number of extrema, and a smaller but still 
vast number of minima. We focus on three specific prob¬ 
lems - the number of minima (relative to the number of 
extrema) that survive when cross-couplings between the 
individual fields are turned on, the average separation 
between extrema in a random landscape, and the infla¬ 
tionary trajectories that exist within the landscape. Cal¬ 
culations with a similar flavor to those presented here can 


^ We use the term “landscape” to refer to any potential with large 
numbers of scalar degrees of freedom - whereas the phrase “string 
landscape” denotes the actual landscape (if any) that is provided 
to us by string theory. 


be found in Arkani-Hamed et al.'s work on the “friendly 
landscape” [10], extended by Distler and Varadarajan to 
cross-coupled fields through the use of random polyno¬ 
mials and modular forms [11]. The latter paper, in par¬ 
ticular, examines a very similar problem to ours, using 
algebro-differential techniques to analyze the minima of 
the landscape. Here, we focus upon the random matri¬ 
ces arising from the Hessians which characterize the ex¬ 
trema of a random cross-coupled potential. We find that 
the cross-coupling terms must be suppressed relative to 
the self-interaction terms by ^ 1/Vn, where N is the 
dimensionality of the landscape, or the number of min¬ 
ima is effectively reduced to zero. This coincides (up to 
a numerical factor) with the separation of scales needed 
to ensure that Newton’s constant is safe from radiative 
corrections [10]. We also look at higher order corrections, 
showing that these can induce an additional log {N) term 
in this result. Further, random matrix theory models of 
the landscape are considered by [12-14], in the context 
of quantum cosmology. 

In addition to estimating the number of vacua permit¬ 
ted by a random multifield potential, we also consider the 
very large number of inflationary trajectories associated 
with a random landscape. Slow roll inflation will typi¬ 
cally occur near a saddle, rather than at a minimum - and 
saddles are much more numerous than minima or max¬ 
ima. Many inflationary models have been built inside of 
string theory (e.g. [15-37]), and a number of these make 
explicit reference to the properties of the landscape. In 
particular, folded inflation [33] and chained inflation [34], 
and rely on the presumed combinatorics of the landscape 
to concatenate several successive periods of inflation and 
we will see that a generic landscape is not likely to have 
the features these models need. Random models of in¬ 
flation have been considered in other contexts. Monte 
Carlo Reconstruction provides a mechanism for generat¬ 
ing an arbitrarily large number of random inflationary 
models, via the Hubble Slow Roll formalism [39]. More¬ 
over, given the extension of the Hubble Slow Roll for¬ 
malism to multifield models, the generalization of Monte 
Carlo reconstruction to more general potentials would 
be straightforward [41] although of limited value given 
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the large number of free parameters this necessarily in¬ 
volves. In separate work, Tegmark has looked closely at 
the distribution of observable parameters associated with 
random inflationary models [40]. Finally, the problem 
tackled here has obvious analogues in condensed matter 
system - see, for example, [42]. 

In Section 2, we review the expected number of ex¬ 
trema in a random landscape, and the proportion of these 
extrema which are actually minima. We then consider 
the impact of cross-couplings on the distribution of min¬ 
ima in two limits. When the cross-couplings are of the 
same order as the self-couplings, we use random matrix 
theory to show that the number of minima is essentially 
zero - or one, if we have a priori knowledge that the 
potential is bounded below. Conversely, for small cross¬ 
couplings we show that if the typical mass-scale in the 
landscape, M, satisfies M < Mpi/^/N the minima are 
safe. In Section 3 we review the inflationary dynamics 
and perturbation spectra of multifield models, and their 
implications for the landscape. We show that if these po¬ 
tentials are bounded in field-space, a weak version of the 
r] problem appears even without reference to supersym¬ 
metry breaking, and that viable models of inflation in a 
landscape are likely to involve either an inflaton rolling 
toward a saddle (hybrid inflation) or several scalar fields 
working cooperatively to drive an inflationary epoch (as¬ 
sisted inflation). In Section 4 we compute the average dis¬ 
tance between extrema in the landscape and show that 
this is typically larger than Mpi, showing even though 
the landscape can contain a vast number of minima, they 
are typically well separated. 


II. EXTREMA OF RANDOM LANDSCAPES 
A. Counting and Classifying Extrema 



FIG. 1: The eigenvalue distribution of a 1000 x 1000 symmet¬ 
ric matrix with entries drawn from a standard normal distri¬ 
bution is shown, with the theoretical distribution plotted for 
reference. 


effectively continuous (assuming they are not massively 
degenerate), yielding a discretuum of values for A.^ In 
this case minima are so closely spaced that no conceiv¬ 
able measurement of A alone could uniquely determine 
which vacuum corresponded to our universe, even if we 
were able to exhaustively catalog all the minima in the 
landscape. 

Recalling elementary multivariable calculus, V pos¬ 
sesses an extremum if all first derivatives vanish at a 
point. Extrema are classified via the Hessian, 


H = 



d^V 

d4>l 
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If H only has positive (negative) eigenvalues, we have 
a local minimum (maximum), whereas if H has mixed 
positive and negative eigenvalues, we have a saddle point. 


Susskind has given a simple geometrical argument that 
the number of minima in the string landscape is most 
likely very large [4]. Consider a landscape with N scalar 
degrees of freedom, 0^, i = 1, - ■ ■, N, and the functional 
form 

m = ( 1 ) 

i 

where the fi each have at least one extremum. Since 
there are no cross terms, any combination of (jji that ex- 
tremizes all the fi yields an extremum of L. If each 
fi has ai extrema, the total number of extrema of L 
is ]([ j ai = a^ , where a is the geometric mean of the 
ai- Continuity requires that every second extremum of 
fi is a minimum, so the landscape has (a/2)^ vacua. 
In the string landscape, N is expected to be on the or¬ 
der of a few hundred. Provided a is not too close to 
unity, ^ 10^^°, tempting one to use anthropic ar¬ 
guments to “explain” why the cosmological constant, A, 
is apparently fine-tuned to this level of precision. More¬ 
over, when 10^^°, the possible values of A are 


B. Extrema and Cross-coupled Fields 

Since E is a function with random coefficients, H is a 
random matrix. The theory of random matrices is rich 
and beautiful, with applications all over physics [45]. The 
simplest case, where the landscape is defined by equa¬ 
tion (1), is trivial: H is diagonal, and the eigenvalues 
are the second derivatives of the fi at a given extremum. 
These are equally likely to be negative or positive, so we 
have 2^ equally weighted permutations of the signs of 
the eigenvalues. Only one of these configurations corre¬ 
sponds to a minimum, and almost any randomly selected 
extremum will be a saddle. Thus, of the extrema, 
/2^ = {a/2)^ are minima, recovering the result of 
the previous subsection. 


^ Although even if A does scan through some range of values, there 
is no guarantee that this range includes zero [10]. 
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FIG. 2: We show the proportion, / of randomly generated 
symmetric N x N matrices for which all eigenvalues have the 
same sign. The curve is well fit by ~ exp(—with the 
key observation being that / is decreasing faster than expo¬ 
nentially with N. 


From a physical perspective, however, there is no rea¬ 
son to expect that the cross-terms in the landscape van¬ 
ish, so let us add simple cross-couplings. Near an ex¬ 
tremum, these take the form 

(3) 

i i^j 

We can now consider the opposite, and equally unphysi¬ 
cal limit, where all the terms in H are of equal size, and 
drawn from identical, independent distributions. In this 
case, we need to turn to a more sophisticated analysis. If 
the Sij are drawn from a normal distribution, our Hessian 
is drawn from the Gaussian Orthogonal Ensemble, and 
the resulting eigenvalues obey the Wigner semi-circle law 
[45]. The density of eigenvalues,^ E{x), is 


ture of this distribution is that it is “rigid” - fluctua¬ 
tions in the eigenvalue distribution are very small. Con¬ 
sequently, the likelihood of a large fluctuation putting all 
the eigenvalues on the same side of zero is vastly smaller 
than 1/2^. 

A great deal is known about the properties of fluctu¬ 
ations away from the semi-circle law (e.g. [46]) and we 
can estimate the likelihood of all the eigenvalues hav¬ 
ing the same sign. Firstly, change variables from Xi to 
= Xi — Xi-i,ai = xi- Without loss of generality we 
can assume that xi < ■ ■ ■ < xn, so the Oi are all positive 
for i > 1 and we can drop the absolute value symbol in 
equation (5). The probability that all the eigenvalues are 
positive is 


p{xi > 0,Vi) 


/“ dai fg°° da 2 --- daNPN{ai) 
/Too ^“1 /o°° da 2 --- /p“ daNPN{a,) 


( 6 ) 

The integrals over ai, i > 2 necessarily have lower bounds 
of zero - whereas all the eigenvalues must be positive if 
ai > 0, and the integral over ai thus starts at —oo in 
the denominator, and 0 in the numerator. Performing 
only the integrals over ai, and leaving the positions of 
the other eigenvalues unspecified, after some algebra and 
a great deal of cancellation we find 


p{a2, - ■ ■ ,aN\ai > 0 ) 


1 

2 


I — erf 



ai+i(iV - 



( 7 ) 

This calculation is exact up to this point. We now guess 
that the average Ui is given by assuming that th e ei gen- 
values are uniformly distributed between (0, V2N) so 
tti ~ y/?^j2. After making this brutal approximation'^ 
and employing the asymptotic form of the error function, 
we find 


E{x) 


^[2N-x^f\\x\<Vm . 4 . 

0, [xi > VT/V 


p{xt > 0,Vi) 


iexp 


(4)f 


( 8 ) 


where the normalization ensures that when E(X) is in¬ 
tegrated over all values of A the total number of eigen¬ 
values is N. Figure 1 shows the eigenvalue distribution 
for a specific 1000^ matrix. For the Gaussian Orthogo¬ 
nal Ensemble, the joint probability distribution for the 
eigenvalues is 

Pn{xi, ■ ■ ■ ,xn) = Givexp)-^ \xj - Xk\ (5) 

j<k 

where Gat is a normalization factor and P gives the likeli¬ 
hood of finding drawing a random matrix with the eigen¬ 
values {xi, • • •, xat} [45]. For our purposes, the key fea- 


In Figure (2), we have plotted p{\i > 0,VI) for ran¬ 
domly generated symmetric N x N matrices. This data 
is accurately fitted by the heuristic relationship 

logp = -0.32491Af"2'°°^®'^. (9) 

Despite the approximation used to obtain equation (8), 
we have confirmed that the likelihood that all the eigen¬ 
values of an A^ X A^ symmetric matrix have the same sign 
scales as e~‘^^ . The measured constant differs slightly 
from —0.25, although given the simplicity of our approx¬ 
imation the agreement is perhaps surprisingly good. 

Tracy and Widom have proved very general results 
about the expected size of the largest eigenvalue, and the 


® These matrices are constructed by taking a random matrix with 
elements drawn from a Gaussian of unit width and then sym¬ 
metrizing. Changing the width of the distribution rescales x. 


^ This is crude for two reasons. Firstly, the ai are not all the same, 
and secondly the maximal eigenvalue is not likely to be exactly 
\/2N after a fluctuation this large. 
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likelihood of fluctuations away from the expected maxi¬ 
mal value [46]. The likelihood of hnding that the largest 
eigenvalue is zero is of course equivalent to the probability 
that all the eigenvalues have the same sign. In principal, 
one could extract an exact formula for this ratio at large 
N from [46], but since we are working with a very simple 
model, a more sophisticated calculation is probably not 
justihed. 

The strong suppression of large huctuations away from 
the semi-circle law can be understood physically. The 
eigenvalues of a symmetric matrix are physically equiva¬ 
lent to the distribution of charged beads along an inhnite 
wire with a quadratic potential centered on the origin, 
and logarithmic interaction potentials between the beads. 
The overall potential causes the beads to cluster near the 
origin, but their mutual repulsion forces them to separate 
[45]. The likelihood of a fluctuation that places all the 
beads on the same side of the origin is very small, and 
vastly less than the 1/2^ likelihood that the eigenvalues 
of a diagonal matrix will all have the same sign. 

The number of extrema in the landscape was calcu¬ 
lated after assuming that the fields are uncoupled, so 
that the number of extrema encountered as one adjusts 
a single held does not depend on the value of the other 
helds. We are now tacitly assuming that even in a cross- 
coupled landscape the number of extrema depends only 
exponentially on N. The precise form of this depen¬ 
dence does not matter, since the chance of any given 
extremum being a minimum (or maximum) now decays 
super-exponentially with increasing N. In this limit, in¬ 
stead of having ^ 10^^° vacua, the potential has none - 
or perhaps just one, if it is bounded below. In a land¬ 
scape with cross-couplings of the same order as the diag¬ 
onal terms, the held point will never be prevented from 
rolling and will naturally how to the global minimum (as¬ 
suming that one exists). While such a potential yields a 
universe that is now in a unique and stable vacuum, one 
cannot use combinatorial arguments to dilute the cosmo¬ 
logical constant problem. 


C. Landscapes With Small Cross-couplings 


coupling terms can be without destroying the expectation 
that 2~^ of the extrema are actually minima. 

We write the Hessian as H = X + U, where X = 
diag(a;i, a: 2 , • ■ •, a^w) is a diagonal matrix with distinct, 
positive eigenvalues xi and U is an N x N symmetric 
matrix.^ The eigenvalues of H, Xi can be written down as 
a perturbation series using Moore-Penrose inverses [47], 
and up to third-order in the elements of U, 


Xi 


Xl + Ull 


7^1 




EE 


UjlUij 

(xj - xi){xj - Xl) 


( 10 ) 


Assume our Hessian is associated with a minimum of 
the diagonal landscape X. If the average eigenvalue of 
X is on the order of 1, then we guess that the lowest 
eigenvalue of X is roughly 1/iV at a local minimum. We 
immediately see that only the lowest lying eigenvalues of 
X are at risk of having their signs changed by adding the 
cross-coupling terms. 

Near a minimum, the Xi are the masses of the N scalar 
degrees of freedom describing oscillations away from that 
minimum. At lowest order, these are given by the xi. The 
diagonal terms of the symmetric contribution, U, corre¬ 
spond to the loop corrections to these masses. To prevent 
the mass corrections from changing the sign of the lowest 
eigenvalue, the standard deviation uu must be (at least) 
several times smaller than 1/N.^ For dehniteness, we ex¬ 
press this bound as ]«] ~ e < l/(3iV) - where the 3 is 
chosen so that a 3cr term in U has the ability to change 
the sign of an eigenvalue."^ 

If we assume the entries of U are induced by Planck 
scale operators, we can relate the dimensionality of the 
landscape to a separation of scales. Restoring units to 
A, A ~ 


m2 
^ m 

M ^ I 
Mpi ~ VW 


( 11 ) 


In reality, it is reasonable to assume that the land¬ 
scape can be decomposed into two pieces - a diagonal 
part corresponding to the /i((/>) which are of order unity, 
and much smaller cross-terms of order e. Physically, the 
diagonal terms are self-interactions and the cross terms 
come from interactions between the individual fields. If 
the characteristic mass-scale of the tree level fi{4>i) is M 
then e has dimensions of mass-squared, and we expect 
e ~ MVM2,, in the absence of special symmetries and 
assuming that the cross-terms are induced by gravita¬ 
tional interactions. 

If e ~ I (that is, the cross-couplings arise from opera¬ 
tors of the same order as those that produce the diagonal 
terms), we have established that almost any extremum 
will be a saddle. We now determine how large the cross 


Interestingly, Arkani-Hamed et al. derived essentially the 
same relationship between the Planck scale and the scale 
of the landscape by requiring that Newton’s constant be 
protected from radiative corrections: they find an ex¬ 
tra numerical factor arising from the constant term in 


® If the eigenvalues of X are not distinct the problem becomes a 
little more complicated, but will not affect out conclusions here. 

® We assume that the entries of U and X are uncorrelated, and 
that all the Uij are drawn from similar distributions. 

^ The “3” here should really be obtained from an expansion of the 
tail of the integrated Gaussian distribution. However, since a 4cr 
fluctuation is 17 times rarer than a 3 <t fluctuation (and the ratio 
increases as we move further from the mean), this is a very mild 
correction and we ignore it. 
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the Einstein-Hilbert action, which we have effectively ab¬ 
sorbed into Mpi. 

At this point, however, we have not examined the con¬ 
vergence of the series in equation (10). Looking at the 
second term, the numerator scales like e^. Assuming that 
the eigenvalues are evenly spaced, so that Xi ~ i/N, we 
estimate the magnitude of the sum as follows: 


S = 

^ S' « 

s ~ 


-E 




{Xi - Xl) ’ 


N 


N 






i=2 N N 


-e^NlogN. 


( 12 ) 

(13) 

(14) 


In order to the protect the lowest-lying eigenvalue, ^ 1/N 
we need 




1 


log N 


(15) 


This term always has a negative coefficient, and tends to 
render the lowest-lying mass term at a given minimum 
negative. We now have a stronger constraint on the terms 
of U than the e ~ 1/N derived above. The off-diagonal 
terms of U now contribute, and in addition to correcting 
the masses they also induce a small rotation in the or¬ 
thogonal modes, which are given by the eigenvectors of 
X + U, not X alone. This result was obtained using an 
explicit ansatz for the eigenvalue spacing of X. Given a 
full understanding of this spacing, which would be possi¬ 
ble in a detailed physical model, it would then be worth 
analyzing these logarithmic corrections carefully. How¬ 
ever, given the generic potentials we are studying here 
we will not pursue this issue further, beyond noting that 
the bounds on e can easily be tighter than 1/A^^. 

The third term of (10) is benign, since the summation 
resembles but the extra factor of N is deleted by 

the uu term in the numerator, and we do not see a log¬ 
arithmic divergence. The final term is again safe, since 
its numerator takes on random signs and its contribu¬ 
tions add incoherently. Moreover, with this assumption 
we see that the above series is at least asymptotically 
convergent. If typical entries in U have magnitude 1/cA^, 
where c < 0{.l), then this series can be safely truncated. 
In Figure 3 we demonstrate these results heuristically, 
showing how preserving the sign of the lowest eigenvalue 
depends on our choice of e. 

We conclude that even a mild separation of scales en¬ 
sures that the off-diagonal terms in the Hessian matrix 
are small enough to preserve the intuition that a land¬ 
scape indeed has multiple minima. We also see that the 
separation of scales needed for the protection of Newton’s 
constant [10] picks up an additional logarithmic term via 
this calculation. For N ~ 1000 this is still less than a full 
order of magnitude, but it is tempting to speculate about 
the connection between this heuristic argument and the 
separation between the GUT scale and the Planck scale. 
The calculation here implies an upper bound, but cannot 
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FIG. 3: We plot the proportion of matrices H = X + U whose 
lowest eigenvalue is negative, where A is a diag(l, ■ ■ ■, N)/N 
and [/ is a symmetric perturbation whose terms have standard 
deviation a. These results are drawn from a sample of 10000 
matrices at each value of a, with N = 100. 


be taken to imply that this bound is saturated in prac¬ 
tice. In particular, there are well-established examples 
(e.g. shift symmetries for axions and other Goldstone 
bosons, and locality in the case of fields associated with 
cycles that are well separated in the compact directions) 
where the e terms are very strongly suppressed. 

The next level of sophistication, which we will not pur¬ 
sue here, would be to analyze a landscape where most of 
the fields were only weakly coupled to each other, but 
with each field coupled to several others at lowest or¬ 
der. This yields a band-diagonal Hessian matrix. As the 
width of the band increased one would move from the 
1/2^ diagonal limit to the exp(—Af^) limit found for the 
fully symmetric case. Given an estimate of the width 
of this band one could then deploy more sophisticated 
random matrix techniques to derive the eigenvalue dis¬ 
tribution for the Hessian. 


III. INFLATIONARY DYNAMICS 

Many discussions of the landscape focus upon the num¬ 
ber of minima, since these are the points that are asso¬ 
ciated with the discrete vacua. Given a vast number of 
minima, we can at least assume the existence of vacua 
where the height of the landscape will match the ob¬ 
served value of cosmological constant, along with any 
other fundamental parameters that are effectively “en¬ 
vironmental.” There is no guarantee that these special 
minima are attractors, but provided they have a non¬ 
zero basin of attraction in the initial conditions space, 
we “explain” the observed value of A with a very mild 
amount of anthropic reasoning. From a cosmological per¬ 
spective, however, we care about both the vacuum itself 
and how the universe arrives there. That is, what is the 
trajectory through the landscape that leads to this final 
(meta)-stable vacuum from a presumably random initial 
configuration? In particular, does this trajectory include 
a region where the height of the landscape changes slowly, 



6 


and we see a sufficiently long period of inflation? In this 
case, the combinatorial arguments advanced to “solve” 
the cosmological constant problem may also address the 
tunings present in most models of inflation. 

Inflation in a landscape (or with any set of scalar fields 
with canonical kinetic terms, for that matter) can occur 
in two different ways: 

• The field point is near a minimum with 1^ > 0, and 
the fields do not possess sufficient kinetic energy to 
climb over the barrier. 


• The field point evolves slowly enough to ensure that 
P < -p/3. 


The first possibility is realized today if the dark energy 
is generated by an actual cosmological constant, rather 
than some quintessence-style scenario where w —1 or 
w ^ 0. The second possibility is the prototype for pri¬ 
mordial inflation, since if the field point is trapped in 
a local minimum it can only escape by tunneling, resur¬ 
recting the bubble problem first recognized in the original 
model of inflation. 

If the Hessian is diagonal, 1 in 2^ extrema will be 
a minimum, but it also follows that there are N in 2-^ 
saddles with a single downhill direction. For successful 
inflation, we need the downhill direction to be relatively 
flat. Near a saddle, the first derivatives necessarily van¬ 
ish, and our criteria for flatness will focus on the second 
derivatives in the downhill directions. 

We will restrict our attention to a flat FRW universe 
~ which is reasonable after some inflation has occurred, 
but avoids questions surrounding the onset of inflation. 
At least locally in the landscape, we can transform the 
fields to ensure that they have canonical kinetic terms, 
so 


d 

a 


SttG 

SttG 


N 

^(</>) + XI 2 

N 

i=l 




dVjcj)) 

d(j)i 


0 . 


(16) 

(17) 

(18) 


Here (() is a shorthand for a vector in fieldspace, (j) = 
{4>i, • • •, 4>n) and labels i run from 1 to N, and H is the 
usual Hubble parameter. 

To sustain inflation, the kinetic term in eq. (16) must 
be much less than the potential term, or 


N 




(19) 


2=1 


We now stipulate 




8irG ^ \ r 

1—1 1—1 \ 





( 21 ) 


by requiring that (j) is small enough to be dropped from 
the equations of motion. Inflation thus occurs near points 
where all the first derivatives are very small or, equiva¬ 
lently, in the vicinity of an extremum.® Consequently, by 
counting the extrema (as opposed to just the minima) 
in the landscape, we can estimate the number of distinct 
inflationary trajectories that can be supported within it. 

In order to produce a phenomenologically acceptable 
period of inflation, the second derivatives of the potential 
must be small, at least in the downhill directions. This 
is ensures that slow roll applies not just at a point in 
fieldspace, but over some finite range, so that the infla¬ 
tionary epoch is not transient. Moreover, the first and 
second derivatives of the fields combine to determine the 
spectral index, of the perturbations - but near a sad¬ 
dle, the spectrum is likely to be dominated by the second 
derivative terms. 

Given our parameter a and the requirement that the 
\(j)i\ < Mpi if the landscape is to be protected from large 
non-perturbative modifications we can guess the average 
value of ry = STr/MpiV"/V, assuming a single evolving 
field. For definiteness, take V ~ sin{2a(j)/TrMpi), 


„ AW 

2.Mpi 4^/2M|,; ■ 

Thus the expected value of r] will be 
ira'^ TT 


( 22 ) 


(23) 


We have taken a>4 in the last step, consistent with the 
requirement that {a/2)^ is a very large number. This 
is a very crude version of the ry problem which plagues 
many supersymmetric models of inflation, although the 
only physical input is the requirement that the field val¬ 
ues remain sub-Planckian. At this level, the rj problem 
does not appear to be especially onerous. A reasonable 
perturbation spectrum needs ry ~ 0.01, and this mild tun¬ 
ing could easily be supplied by relying on the vast num¬ 
ber of potential inflationary paths.® In realistic models 
of inflation, ry is derived from calculable supersymmetric 
corrections to the potential, and these must be analysed 
more carefully. 

There is a hidden assumption in the analysis above: 
that we are only moving in a single direction. However, 
if a landscape extends ~ Mpi in any given direction, 


® Strictly, if the first derivatives are all very small near a point, 
this does not guarantee there is an extremum nearby. However, 
this qualificiation does not change the essence of our argument. 

® One of us [RE] will examine the likely values of rj in single field 
inflation in a forthcoming paper [50]. 
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Pythagoras’ theorem tells us that it measures ~ \/NMpi 
from corner to corner. Consequently, if the structure of 
the landscape permits the field point to roll diagonally, it 
may traverse a distance far greater than Mpi. This is pre¬ 
cisely the scenario envisaged by assisted inflation, where 
multiple fields act coherently to emulate a single min¬ 
imally coupled scalar field [9]. This effectively embeds 
inflationary models with a single scalar field and trans- 
Planckian vev within the landscape. However, in order 
to see this, we must transform to a new basis in which 
the only downhill degree of freedom runs from one corner 
of the landscape to another. If we view the landscape as 
a purely random potential, assisted inflation does not ap¬ 
pear natural (the position taken in [33]), since it requires 
a correlation between the masses of a large fraction of 
the fields, and this reduces the effective dimensionality 
of the landscape, undermining the strength of combina¬ 
torial arguments. In particular N-flation [38] proposes a 
mechanism for realizing this type of potential within the 
landscape. 


IV. DISTANCE BETWEEN EXTREMA IN THE 
LANDSCAPE 

We now ask whether the set of extrema is large enough 
so that any random point will be “near” a minimum or 
saddle, and the minima and extrema thus form a dense 
set inside the landscape. In the converse case, extrema - 
while numerous - still correspond to very special points 
inside the landscape. We make considerable progress by 
appealing to the Central Limit Theorem (see e.g. [49]). 
We assume the (pi at each extremum are drawn inde¬ 
pendently from identical, flat distributions but we could 
relax this assumption without undermining our results. 

Take an arbitrary reference point in the landscape, 
(j) = (^ 1 , • • •, (Pn), and consider a set of M randomly dis¬ 
tributed points in the landscape (p^ = {p\, • • •, where 
j G {I,---,M}. In practice these points will be the 
set of minima or inflationary saddles. Assume a simple 
Euclidean metric on the landscape, so that the distance 
between p and an arbitrary obeys 

N 

i=l 

The Central Limit Theorem tells us that X, a sum over 
N random variables Xi each with mean fii and variance 
CTi tends to a Gaussian as N becomes large. Constructing 
X as follows 


A = 


N N 


E N 2 


(25) 


gives a Gaussian with unit variance and zero mean. Con¬ 
sequently, if the Pi are randomly distributed, the distri¬ 
bution of approaches a Gaussian random variable. 


We draw the pj from a uniform distribution of width 
To keep the field values bounded by Mpi, we need 
A ^ Mpi. The mean and variance of each term in the 
sum in equation (24) is, 

^ ^ dpl[{pl - pif - y.iY 

= 1^2^2^1a4. (27) 

Putting this together, 

+ ■»! , 28 ) 

where is a normalized version of the distribution of 
. 

The above expression applies to a specified reference 
point in the landscape, p. If we also choose p at random 
from the same distribution as the pP we have the dis¬ 
tribution of distances between any two randomly chosen 
points. The means and variances in equations (26) and 
(27) become two dimensional integrals over both pj and 
pi, yielding 


(r^) = (29) 

cr(r2) = A^y^V (30) 

In the large N limit, the distribution for the square of 
the distance between two randomly chosen points is 


P(r^) 


1 

/ 45 

45(r2 - 2f A2)2 

A2\ 

' 567rV ^ 

56AA4 


(31) 


In terms of , a separation of A^ corresponds t o a set of 
points that lies 2VA^/3 — A^ « 2VA^/3 or ^/5NJ7 « 
a/TV standard deviations from the mean. In fact, this 
far out into the tail we cannot trust the Gentral Limit 
Theorem result: the tails of P(r^) are suppressed relative 
to the Gaussian form, since 0 < < fVA^ while P(r^) 

computed from (31) is non-zero outside these regions. 

Instead, we can estimate P(r^) geometrically when 
is far below the mean. A hypersphere of radius S sur¬ 
rounding a point has volume 

27r^/2 


Only the mean and variance of the distribution will enter the 
calculation, so choosing a different distribution will not change 
the results here dramatically. 



where we are allowing the volume of interest to have a 
radius S which differs from the size of the box that con¬ 
tains the landscape. Conversely, the landscape itself is 
(by assumption) a cube with volume 

Cn = (2A)^ . (33) 

Taking the ratio of these two quantities, the odds of ran¬ 
domly choosing two points in the landscape and find¬ 
ing that they are separated by a distance is 

Sn/Cn- This quantity rapidly approaches zero as N be¬ 
comes large for fixed A, thanks to the familiar result that 
the N-cube is vastly larger than the N-sphere. Quantita¬ 
tively 


P(r2 < 6 ^) 


Sn 

Cn 


1 

V ttN 



(34) 


Moreover, this is actually an overestimate, since it as¬ 
sumes the whole A-sphere fits inside the A-cube. 

Recalling that we have (a/2)^ minima, we can there¬ 
fore form (aj2)'^^ j2 distinct pairs of minima. In this 
case, the expected number of pairs of minima whose sep¬ 
arations are less than the characteristic scale of the land¬ 
scape, A is 


A/'(r^ < A^) 


V2^ [2N 



(35) 


Thus we can quickly deduce: 

a<2Ai/4 ^ A/'(r2 < A^) < 1. (36) 


If the inequality in (36) is satisfied, all the vacua in the 
landscape will be separated by a distance greater than 
A. If we had looked for genuinely close points (since A is 
a substantial distance within the the landscape) the con¬ 
straint would have been correspondingly tighter. Numer¬ 
ically, if A ~ 100, a>6 will give a dense set of minima, 
whereas with A ~ 500 would need a>9. Consequently, 
it seems likely that we can regard all the extrema as be¬ 
ing well separated from one another, and thus forming a 
very sparse set within the landscape. 

We can also ask about the likely separation between 
inflationary saddles. In general, there saddles in the 
landscape, but if we want to concatenate two periods of 
inflation without a significant interruption, we want to 
end one trajectory and find that another saddle is within 
a distance H inside field space. Since we have argued that 
the energy scale of the landscape is necessarily reduced 
by a factor of 1/VN relative to the Planck scale, H ~ 
5 ~ M'^/Mpi, and S/A // M’^/Mp^ ~ 1/A. Inserting 
these values into equation (34) 


P(r2 < p2) = 


VttA [ 


e I 
2N1^ 


-I JV/2 


(37) 


and the likelihood of connecting two inflationary sad¬ 
dles by chance drops super-exponentially as A increases. 


Even though the number of extrema grows exponentially 
with A, the ratio volumes between the hypercube and hy¬ 
persphere grows super-exponentially with A. Moreover, 
the separation of scales needed to protect the minima 
grows with A, and since this is inversely proportional to 
the scale of inflation, the likelihood of finding two adja¬ 
cent inflationary trajectories is further suppressed. This 
conclusion applies to both “chained inflation” [34] and 
“folded inflation” [33]. Since we are looking at a purely 
random potential, it is possible that both of these models 
might be realized in the context of the actual string land¬ 
scape. On the other hand, multistage models of the form 
discussed in [35] are safe, since the inflationary stages are 
well separated. However, they are not rendered natural 
simply by the large number of extrema that the land¬ 
scape contains, but must be justified by showing that in 
the actual string landscape there are well-defined sub¬ 
spaces where the density of extrema is much higher than 
in the rest of the volume. This problem is particularly 
acute for chained inflation, as the successive tunneling 
events it needs will be strongly suppressed if the minima 
are well separated as the barrier between them is then 
likely to be high. 


V. CONCLUSIONS 

The purpose of this paper has been to explore what 
follows solely from the large dimensionality of the string 
landscape, rather than from its detailed mathematical 
properties. 

As was stressed earlier, the results here apply to gen¬ 
eral “random” potentials, but it is entirely possible that 
the actual distribution of minima in the string landscape 
is sufficiently non-random to thwart the arguments made 
in this paper. However, these results will let us distin¬ 
guish between intrinsically stringy effects, and phenom¬ 
ena expected to arise in any physical model with a large 
number of scalar fields. Even at this point, we can iden¬ 
tify a number of ways in which string theory may yield a 
distribution of vacua which differs from those posited by 
our analysis. For instance, [51] examines no-scale mod¬ 
els where the potential is non-negative and the Hessian 
matrices positive semi-definite. There are flat directions 
at tree-level but no tachyons, leading to a set of Hessian 
matrices that differs sharply from those considered here. 
After SUSY breaking, the extrema seen in these mod¬ 
els [7] still have Hessian matrices whose detailed struc¬ 
ture ensures that the fraction of extrema which are ac¬ 
tually minima exceeds 1/2^. Finally, we can seek out 
sub-regions within the landscape which contain an extra 
scale M*, so the separation of vacua in field space is con¬ 
trolled by M* rather than Mpi. Conversely, we might 
well find that variants of the arguments developed here 
still apply within the subspaces of the landscape in which 
the vacua are “over dense” or otherwise atypical. 

Interestingly, while the number of minima produced by 
the simplest model - a large set of non-interacting fields 
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- can be estimated trivially, turning on cross-coupings 
quickly changes the picture. In particular, if the cross 
terms between the fields have the same size as the diago¬ 
nal quadratic terms, we can use random matrix theory to 
show that almost all extrema of the landscape are saddle 
points, and that there will be few or even no minima. In 
the limit where the cross-terms between the fields are nat¬ 
urally suppressed relative to the self-interaction terms, 
we can quantify the amount of suppression required to 
ensure that the minima are not removed by these inter¬ 
actions. 

Finally, the distances between minima in the landscape 
are typically substantial, even though the number of min¬ 
ima grows exponentially with the dimensionality of the 
landscape. This has immediate consequences for mod¬ 
els such as chained inflation [34] or folded inflation [33], 
which concatenate several separate periods of inflation. 
These models are not ruled out, but the necessary com¬ 
binatorics cannot be provided simply by an appeal to the 
dimensionality of the landscape. 

We see several possible extensions of the approaches 
outlined in this paper, and feel that the appearance of 
random matrix theory in this context is particularly sug¬ 
gestive. In particular, one of the authors [RE] and Liam 


McAllister are currently using random matrix techniques 
to analyse the mass-spectrum found in models of N- 
flation [38]. The extent to which the properties of the 
actual string landscape matches the expectations derived 
solely from the large-N considerations explored here re¬ 
mains an open question. However, we have shown the di¬ 
mensionality of the landscape alone is sufficient to draw 
conclusions about the likely set of extrema that exist 
within it. 
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